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ABSTRACT 

We propose to use spatial correlations of the kinetic Sunyaev-Zeldovich (KSZ) flux as 
an estimator of the peculiar velocity power spectrum. In contrast with conventional 
techniques, our new method does not require measurements of the thermal SZ signal 
or the X-ray temperature. Moreover, this method has the special advantage that the 
expected systematic errors are always sub-dominant to statistical errors on all scales 
and redshifts of interest. We show that future large sky coverage KSZ surveys may 
allow a peculiar velocity power spectrum estimates of an accuracy reaching ~ 10%. 



Subject headings: cosmology: large scale structure; theory- 
cosmic microwave background; cosmology: observations; 
cosmology: theory; dark matter; distance scale 



1 INTRODUCTION 

Peculiar velocity-distance galaxy surveys have provided us 
with maps of line-of-sight peculiar velocity fields, w(x) = 
V ■ x./x, where V is the peculiar velocity at position x. 
Such surveys were successfully used in the past to con- 
strain power spectra of initial density fluctuations and to 
estimate the cosmological density parameter, Q m (Vitto- 
rio, Juszkiewicz & Davis 1986; Groth, Juszkiewicz & Os- 
triker 1989; da Costa et al. 1998; Feldman, et al. 2003a; 
Feldman, et al. 2003b; Sarkar, Feldman & Watkins 2007; 
Watkins & Feldman 2007; Feldman & Watkins 2008). For 
an overview of the subject, see (Strauss & Willick 1995) 
and (Courteau & Willick 2000). However, all these velocity- 
distance surveys are limited by errors in the range of 15 to 
20% of the estimated distance. Measurements based on type 
la supernovae have slightly lower error per source, due to 
smaller dispersions in type la supernovae intrinsic luminosi- 
ties (see, e.g. (Hui & Greene 2006; Gordon, Land & Slosar 
2007) and references therein for theory aspects and (Bon- 
vin, Durrer & Kunz, 2006; Haugboelle et al. 2007; 1; Wang 
2007; Watkins & Feldman 2007) and references therein for 
recent measurements). In contrast, for cluster line-of-sight 
peculiar velocities, v, derived from the kinetic Sunyaev Zel- 
dovich (KSZ) effect, the errors grow much less rapidly with 
distance. KSZ cluster surveys may open new possibilities for 
studying large-scale flows (Haehnelt & Tegmark 1996; Kash- 
linsky & Atrio-Barandela 2000; Aghanim, Gorski & Puget 
2001; Atrio-Barandela, Kashlinsky & Mucket, 2004; Holder 



2004) and cosmology (Bhattacharya & Kosowsky, 2007; Fos- 
alba & Dore 2007). The currently established method of re- 
covering v from the SZ data requires extra measurements 
of the cluster temperature and its Thomson optical depth, 
t. Other complications include the relativistic thermal Sun- 
yaev Zeldovich (TSZ) effect and the KSZ signal generated by 
internal motions, which is added to the signal from cluster's 
bulk velocity, resulting in a limit in accuracy in determining 
v, of the order of ~ 200 km/s (Knox, Holder & Church 2004; 
Aghanim, Hansen & Lagache, 2005; Diaferio et al. 2005). 

In this paper we propose to estimate the spatial two- 
point correlation function of the KSZ flux from the data and 
then infer the velocity power spectrum. Contaminations to 
the cluster KSZ flux have various clustering properties that 
can be applied to disentangle the signal from its contam- 
inants (see also (Zhang et al. 2004)). We will show that 
at z}t 0.3, the systematics become sub-dominant and the 
statistical errors, at ~ 10% level for South Pole telescope 
(SPT 1 ) (Ruhl et al. 2004), dominate. Throughout this let- 
ter, we adopt the cosmology with Q m = 0.3, £Ia = 1 — fim 
and a normalization parameter erg = 0.9. 



2 THE FLUX POWER SPECTRUM 

The KSZ cluster surveys directly measure the sum of the 
cluster KSZ flux Sksz and various contaminants, such as 
intracluster gas internal flow, radio and IR point sources, 
primary cosmic microwave background (CMB), cosmic in- 
frared background (CIB), etc. The signal is (Sunyaev & Zel- 
dovich 1980) 

1 http:/ /spt. uchicago.edu/, specifically, see the SPT White Paper 
at this site. 
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_ dB v (T) 

>3KSZ — qj, ^KSZ 



(1) 



where B V (T) is the Planck function. At v ~ 217Ghz, 
dB/dT = 540 Jy sr" 1 fiK^ 1 . The KSZ temperature fluc- 
tuation can be expressed as 



ATksz = 2.7K (r) (v/c) = 9fiK v 100 M/0.01, 



(2) 



where (r) is the optical depth averaged over the solid an- 
gle ui, while «ioo = w/(100km/s). We will also use S'ioo = 
<S , KSz(«)/<S , KSz(100km/s). These quantities are, respectively, 
the "normalized" peculiar velocity and KSZ flux, scaled to 
v — 100km/s). The value (r) = 0.01 is typical for clus- 
ter models considered in numerical simulations (see e.g. 
(Yoshida, Sheth & Diaferio 2001) and references therein). 
Since at v ~ 217 Ghz, the non-relativistic part of the TSZ 
effect, which is one of the major contaminants of the KSZ, 
vanishes (Zel'dovich & Sunyaev 1969), we focus on this fre- 
quency. Optical follow-up of KSZ surveys such as dark en- 
ergy survey 2 will measure cluster redshift z with uncertainty 
J; 0.005. 3 The z information allows the measurement of the 
3D auto correlation of measured cluster KSZ flux 



£s(r) = (ASiASj) = (Sksz^Skszj) + ■ 



(3) 



where ASi — Si — IcOi , Si and uji are the total flux and 
the solid angle subtended by the i-th cluster; r = |x< — 
Xj| is the comoving separation of a pair of clusters, and 
/ = ^2 Si I ^2 uii is the mean flux intensity. The total flux 
is a sum of the KSZ signal, Sksz.i, and various noise terms, 
which we study below. 

The correlation (Sksz.iSkszj) is cluster number den- 
sity weighted and thus 

(Sksz.^ksz,,) « ((1 + 5 c {Mi))vi (1 + Sc{Mj))vj) , (4) 

where 5c (M) = b n (M)5 is the cluster number overdensity, 
while 8 is the matter overdensity. Here we use a simple bi- 
asing model, where the parameter b n is determined by the 
cluster mass M (Mo & White 1996) (see also refs. (Kash- 
linsky 1987; Kashlinsky 1991)). The quantity above can be 
written as 



((1 + 8 i )v l (1 + Sj)vj) = (viVj) + (SiSjViVj) , 



(5) 



where the third order moments cancel because of statisti- 
cal isotropy and homogeneity which reflects the symmetry 
expected in the large scale pattern of the cluster distribu- 
tion: (viVjSi) = —(viVjSj). The fourth moment can always 
be written as a sum of products of second moments of all 
possible pairs plus the irreducible or connected fourth mo- 
ment: (1234) = (12}(34) + (13) (24) + (14) (23) + (1234} c . 

The signal we are interested in, (5ksz,«S'ksz,j), is di- 
rectly related to the peculiar velocity logarithmic power 
spectrum A 2 . Here and below, Aq = Pq (fc, z)k 3 /(2n 2 ), 
where k is the comoving wavenumber, z is the redshift and 
Pq is the proper power spectrum of the random field Q. The 
variance is (Q 2 ) = J Q dkAg/k. 

The nonlinear scale at z ~ 1 is k ~ 0.5fa/Mpc, hence for 
wavenumbers k < 1/i/Mpc we can apply the weakly nonlin- 
ear perturbation theory. We also assume here that to first 



2 http://www.darkenergysurvey.org/ 

3 The photo-z of each galaxy has dispersion ~ 0.05. Clusters have 



order, the density and velocity are Gaussian random fields, 
related by the expression 



A v (k,z) = f HaA m / H kVs, 



(6) 



where a is the scale factor, / = d In D/d In a, H = H(z) 
is the Hubble parameter at redshift z — l/o — 1 and 
Ho = H(0), and D(z) is the usual linear growth factor for 
density perturbations. A^ is the logarithmic power spec- 
trum of the dark matter density fluctuations derived using 
the standard transfer function (Bardeen et al. 1986) and 
assuming f2 m = 0.3, Qa = 0.7, a s = 0.9 and Ho = 70 
km/s/Mpc. Under these assumptions, to lowest order, the 
connected fourth moment vanishes because of Gaussianity 
while the nonlinear, mode coupling contribution to the sig- 
nal is well approximated by the expression (Ma & Fry 2002) 



AL(k) « A 2 m (k,z)a 2 v (z) /2 



(7) 



where a 2 — J °° dkA 2 /k is the line-of-sight velocity disper- 
sion. We also assume that the two lines of sight are almost 
parallel. The accuracy of this assumption is better than 
~ 1% for angular separations < 10°. 

Bulk flows are induced by the large scale gravitational 
potential, so the velocity field decouples from the small scale 
density field and the last expression in Eq.(7) remains valid 
and agrees with N-body simulations in the nonlinear regime, 
at redshifts Z < 1 and wavenumbers k > 1/i/Mpc if the 
linear A^ is replaced by its nonlinear analog (Ma & Fry 
2002). 

To express these results in terms of the cluster number- 
density fluctuations 8c (M), we will again turn to a simple 
Press-Schechter prescription, described earlier. The KSZ sig- 
nal is then 



ALz(fc)=4(A?(fc)+6„ 2 A 2 ,) , 



(8) 



> 



100 galaxies and thus the determined z dispersion is J$ 0.005 



where 

p oo 

s N = b"(M)S(M)(dn/dM)dM, N = 0,1, 

J M, 

and b n = si/so- The quantities A„ and A^ are the loga- 
rithmic power spectra of spatial correlation functions {V1V2) 
and (S1V1V262}, introduced above. To describe the density 
field, we use the cluster mass function, written as dn/dM, 
and M- t is the lower mass threshold of SZ selected clusters. 
For SPT, M* ~ 10 14 M Q . The redshift dependence of M» 
(Holder et al. 2000) is weak and we assume a constant M*. 
There is no biasing assumption for velocities because we as- 
sume that the galaxies trace the dark matter velocity field, 
in agreement with observations (Feldman, et al. 2003a). To 
calculate the effective cluster bias, we follow the halo model 
of (Mo & White 1996): b n (M) = l + (v 2 - 1)/(S C /D), where 
5 C ^ 1.686, v = 8 c /(Da(M)) and D(z) is the linear density 
growth factor. The quantity a(M) is the rms density fluc- 
tuation in a spherical volume containing mass M and we 
have ct(M») ~ as = 0.9. For the present estimate we use the 
approximation 6„ = &„(M„). At present time &„(M*,0) w 3. 
This value is consistent with most recent measurements of 
the correlation function of optically selected galaxy clusters 
in the Sloan Digital Sky Survey (Estrada, Sefusati & Frie- 
man 2008). This is, indeed, a simplified version of the biasing 
scheme, however for the scales of interest in this paper it is 
justified. 
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Figure 1. The normalized KSZ cluster flux power spectra at 
z = (upper black solid line) and z = 1 (lower blue solid line) re- 
spectively. The power spectrum, A^.g z (fc) is expressed in dimen- 
sionless units, defined in §2. The dashed and dotted lines show the 
contributions from the vv and vSv8 terms, respectively. The solid 
line shows the sum of these two terms, cr 2 = o" 2 /(100 km/s) 2 . 



The normalization is A? KSZ /s^ 



In Fig. 1 we show the normalized KSZ flux power spec- 
trum at z — and at z = 1. It is clear that the vv contri- 
bution dominates at small k whereas the SvSv term domi- 
nates in the high k tail. A 2 peaks at k ~ 0.05ft/Mpc, be- 
cause of the k 2 denominator, while A^, associated with con- 
taminants traced by point sources, keeps decreasing toward 
large scales. On the other hand, A^sz increases with k be- 
cause of the presence of the mode-coupling term. Note also, 



A KS z has a weak z dependence at 



that compared to A„ 
k< 0.03/i/Mpc. 

We have not considered the effect of redshift distor- 
tion. It reduces the power of the A 2 part by a factor of 
~ exp[— k 2 a 2 /H 2 {a)\. Since a v ~ 200 km/s, the redshift 
distortion can be neglected at k^S 0.2/i/Mpc. On the other 
hand, the redshift distortion increases the power of the A 2 e 
part in the linear regime (Kaiser 1987). Since A 2 S dominates 
over A 2 at 0.1/i/Mpc (see Fig. 1), neglecting the redshift 
distortion tends to underestimate the signal A K sz and thus 
overestimates the relative errors. 



3 NOISE SOURCES 

To study the detectability of A| sz , let us consider the South 
Pole telescope experiment 6 , planned for the near future. The 
SPT will cover the 217 Ghz range and about 4000 deg 2 with 
arc-minute resolution. It is expected to observe more than 
20,000 clusters with masses greater than 10 14 M Q for the 
thermal SZ signal. Virtually all the TSZ clusters have flux 
measurements at 220 Ghz band, a fraction of this flux comes 



from KSZ effects. The corresponding KSZ flux may be too 
weak to detect cluster by cluster, however, our method pro- 
vides a way to extract this KSZ signal at 220 Ghz band by 
statistically weighting all fluxes. Thus, cluster finding is done 
by combining all bands (90, 150, 220 and 270 Ghz). Given 
these TSZ selected clusters, we can measure the CMB flux 
enclosed in the same region at 220 Ghz. The measured flux is 
the sum of KSZ, residual TSZ, CMB and various other con- 
taminations. Indeed, for many individual clusters, the KSZ 
signal is overwhelmed by other contaminants, however, the 
correlation approach described here is able to separate the 
KSZ signal from various sources of noise. 

In the following we discuss the three dominant error 
sources in the following subsections. We assume redshift bin 
size Az = 0.2 and k bin size Afc = OAk. The statistical 
errors scale as Afe -1 ^ 2 . 



3.1 Diffuse foregrounds and backgrounds 

Clusters have typical bulk velocity around several hun- 
dred kilometer per second, so the typical KSZ signal is 
ATksz ~ 20/Ltir(r)/0.01. Primary CMB, which is indistin- 
guishable from the KSZ signal in frequency space, has in- 
trinsic temperature fluctuation AT ~ 100^tK. The cosmic 
infrared background (CIB) has mean T ~ 20/^K, if scaled 
to 217 Ghz (Fixsen et al. 1998). A cluster KSZ filter can be 
applied to remove the mean backgrounds while keeping the 
KSZ signal. This filter must strongly match the cluster KSZ 
profile while having zero integrated area. Since the angular 
sizes of clusters are usually close to several arc-minutes, such 
a filter peaks at multipole I around several thousands and 
thus removes the dominant CMB signal, concentrated near 
l^S 1500. For such filters, at ~ 217 Ghz, the galactic syn- 
chrotron, dust emission, free-free foregrounds, the radio and 
TSZ backgrounds are all negligible due to their frequency or 
scale dependence (see, e.g. (Wright 1998; ?)). Thus we only 
discuss the contaminations of the primary CMB, CIB and 
background KSZ. 

The optimal filter can be constructed from the intr- 
acluster gas profile inferred from the TSZ survey. With- 
out loss of generality, we choose the electron density pro- 
file as n e (r) oc (1 + r 2 /r 2 )" 1 and a matched filter We = 
6(£/£f) 2 exp( — (£/£f) 2 ). For these particular choices, the fil- 
tered KSZ temperature, ATksz, peaks at if « 1.1/9 C w 
3800(1 /6 C ), where 8 C is the angular core radius, and the 
peak value is ATksz ^ ATksz ^ 9v W o[{T}/0M]fiK. We 
adopt this If and 8 C = 0.4/i -1 Mpc/x(z) to estimate the 
noise. Here, x( z ) i s t ne comoving angular distance. 

The correlations of filtered backgrounds (with zero 
mean flux), originating from both their intrinsic and cluster 
correlations, are 



i„(r) « [fu;(M)dn(M)] : 



l + bn ^m(r) 



M0) 



(9) 



where 9 — r/x, while Wb(6) is the corresponding (filtered) 
background angular correlation function and u(M) is the 
solid angle subtended by a cluster of mass M at redshift z. 
Finally, £ m is the dark matter correlation function. 

The same cluster survey measures both J Ludn and 

b n 2 £ m , while the CMB Ce's are known from observations 
(Hinshaw et al. 2003). Alternatively, for the cosmological 
model assumed here, the expected Ct's can be also calcu- 
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Figure 2. The dominant errors of the cluster KSZ power spec- 
trum measurement. The solid lines are the errors caused by the 
cosmic variance and shot noise. The short dash lines are the sys- 
tematic errors caused by the CIB and KSZ background. The long 
dash lines are the statistical errors caused by the primary CMB. 
The dot lines (which are visible only in the upper left panel) are 
the statistical errors caused by sources associated clusters. Sys- 
tematic errors are sub-dominant at fcj; lh/Mpc. 



resolved source subtraction, one is likely able to subtract 
much of these contaminations. Rather conservatively, we as- 
sume that, at v ~ 217Ghz, the total flux contributed by the 
sources associated with clusters is less than ~5x 10 3 Jy/sr 
at z — 0.5 and scale it to other redshifts assuming no in- 
trinsic luminosity evolution. At 217Ghz, the non-relativistic 
TSZ vanishes. The relativistic correction of cluster TSZ ef- 
fect shifts the cross over point slightly to higher frequency 
and thus in principle introduces a residual TSZ signal in 
~ 217Ghz band. But the correction to AT is generally only 
several percent of the TSZ at Rayleigh-Jeans regime (Steb- 
bins 1997; ?). 

The mean flux of the cluster sources, S c \, is subtracted 
in our estimator (Eq. 3). Since the cluster thermal energy, 
IR and radio flux should be mainly determined by local pro- 
cesses, one can omit the possible large-scale correlations of 
these quantities. Thus these sources do not cause systematic 
errors. But since the mean square of the flux fluctuations, 
caused by clusters, {SS^}, does not vanish, they do con- 
tribute a shot noise, which will be estimated in §3.3. Fur- 
ther, they contribute a statistical noise term proportional 
to the cluster clustering, whose expected fractional error is 
T)(k) « 2nA k {SSci)/Sf 00 B k VVk^Ak , for a survey of vol- 
ume V. This can be rewritten as 



5 x 10" 



ZJ 100 



(SSti) AzAk/k 



0.2 0.4 



(11) 



It is reasonable to assume that (SSl) 1/2 < (S c i). Thus the 
error caused by the sources associated with clusters is neg- 
ligible at all scales and redshifts (See Fig. 2). 



lated from first principles (Seljak & Zaldarriaga 1996). Thus 
the CMB contamination can be easily subtracted from the 
correlation estimator, leaving the statistical errors from in- 
trinsic CMB fluctuations over cluster regions as the only 
CMB contaminant. 

In principle, the CIB and KSZ contaminations can 
be removed as well. But since both the amplitude and 
shape of CIB and KSZ power spectra are highly uncer- 
tain, we do not attempt to subtract their contribution 
from the correlation estimator. The CIB power spectrum 
is cf ] i 2 /{2%) ~ (4^tK) 2 (£/10 3 ) ' 7 and the KSZ power spec- 
trum is Cf ) £ 2 /(2ty) ~ (2.7^iK) 2 (Knox, Cooray & Eisenstein 
2001). The resulting upper limit on the fractional systematic 
error is 



(10) 



27rB fc [9((r)/0.01)/iif] 2 ' 
where J = 1, 2 is the summation index; the spherical har- 
monic number, I = k\, is fixed, At = b n A^, while 
Bk = A 2 100 + b n A 2 4100 . Where A 2 ,s 100 is the A 2 , term, 
expressed in units of (100km/s) 4 . Clearly, rj(k) is at most 
several percent at k^, 1/i/Mpc (see Fig. 2). 

3.2 Sources associated with clusters 

The map filter does not remove the contaminations from 
sources associated with clusters. The flux of cluster radio 
and IR sources is ~ 10 3 Jy/sr at 217 Ghz (Aghanim, Hansen 
& Lagache, 2005). Using multi-frequency information and 



3.3 Cosmic variance and shot noise 

The signal intrinsic cosmic variance dominates at large 
scales. The number of clusters is limited, so the shot noise 
is large, even in the linear regime. We believe that there are 
five dominant sources of shot noise. The shot noise power 
spectrum is then given by 



ishot — 



2fnv 2 



E 



-2 



(12) 



where the summation index J corresponds to mean square 
noise values of discrete sources, associated with clusters 
(J = 1), the instrumental noise, the residual flux fluctu- 
ations of primary CMB, CIB, and the background KSZ 
(J = 2 to 5, respectively), n is the mean number density 
of observed clusters, which can be calculated given the halo 
mass function, the survey specification and the gas model. 
We use the expression n(z) = 3x 10~ 5 /(l+2) 3 (/i/Mpc) 3 . We 
estimate the contributions of the five noise sources discussed 
here as follows: b\ w 5[iK and dj w 20/j.K for J = 2 — 5. The 
filtered a of intracluster gas internal flow (Nagai, Kravtsov 
& Kosowsky 2003) is < 10fiK, due to its random nature. 
Lack of knowledge on its spatial distribution, we neglect 
this term. For SPT, the errors caused by the cosmic vari- 
ance and shot noise (~ 10%) dominate over all other errors. 
The systematic errors are virtually always sub-dominant. In 
this sense, our method has the special advantage to measure 
the cluster peculiar velocity. For a future all sky survey, total 
error can be reduced to several percent level. 
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4 CONCLUSIONS 

We presented a new method to estimate the peculiar veloc- 
ity power spectrum using the KSZ effect. The signal (cluster 
bulk velocity correlation) is amplified by direct auto correla- 
tion of the cluster KSZ flux. Consequently, many systemat- 
ics, such as internal flow involved in the usual v p inversion, 
vanish and the majority of remaining systematics can be 
simply subtracted. The correlation method has the special 
advantage that statistical error dominates over systematics 
at effectively all k and z range. 

Ideas similar to those expressed here about prospects 
of probing the peculiar velocity field of clusters at high red- 
shift through KSZ flux correlations have been recently pre- 
sented in (Hernandez-Monteagudo et al. 2006). Our papers 
are complementary: we focus on spatial correlations, while 
(Hernandez-Monteagudo et al. 2006) study angular correla- 
tions. 

Our calculations show that the SPT can measure 
to ~ 10% accuracy at k ~ 0.02-0.3h/Mpc and a future all 
sky survey can improve this measurement by a factor of few. 
These results should be regarded as preliminary - the goal of 
this paper was to present a new idea and a rough estimate of 
the expected signal-to-noise ratio. A comparison with more 
realistic numerical simulations should be the next step. 



5 ACKNOWLEDGMENTS 

PJZ thanks the support of the DOE and the NASA grant 
NAG 5-10842 at Fermilab at the early stage of this work and 
the National Science Foundation of China grants 10673022 
and 10533030. HAF was supported in part by the University 
of Kansas GRF and by a grant from the Research Corpo- 
ration. RJ was supported by (to be added just before sub- 
mission) . AS is supported by the DOE and the NASA grant 
NAG 5-10842 at Fermilab. 



Fixsen, D. J., Dwek, E., Mather, J. C, Bennett, C. L., Shafer, R. 

A., 1998, ApJ, 508, 123 
Fosalba, P. & Dore, O., 2007 PRD, 76 103523 
Gordon, C, Land, K. and Slosar, A., 2007, PRL, 99, 081301 
Groth, E.J., Juszkiewicz, R., & Ostriker, J. P., 1989, ApJ, 346, 

558 

Haehnelt, M.G., & Tegmark, M., 1996, MNRAS, 279, 545 
Haugboelle, T., Hannestad, S., Thomsen, B., Fynbo, J., Soller- 

man, J., & Jha, S., 2007, ApJ, 661, 650 
Hernandez-Monteagudo, C., Verde, L., Jimenez, R. & Spergel, D. 

N. 2006, ApJ, 643, 598 
Hinshaw, G., et al. 2003, ApJS, 148, 135 

Holder, G. P., Mohr, J. J., Carlstrom, J. E., Evrard, A. E., Leitch, 

E. M., 2000, ApJ, 544, 629 
Holder, G., 2004, ApJ, 602, 18 
Hui, L., & Greene, P. B., 2006, PRD, 73, 123526 
Itoh, N., Kohyama, Y., & Nozawa, S. 1998, ApJ, 502, 7 
Kaiser, N., 1987, MNRAS, 227, 1 
Kashlinsky, A., 1987,ApJ, 317, 19 
Kashlinsky, A., 1991, ApJ, 376, L5 

Kashlinsky, A., & Atrio-Barandcla, F., 2000, ApJ, 536, L67 
Knox, L., Cooray, A. & Eisenstein, D., 2001, ApJ, 550, 7 
Knox, L., Holder, G., & Church, S., 2004, ApJ, 612, 96 
Ma, C.-P., & Fry, J.N., 2002, Phys. Rev. Lett., 88, 211301 
Mo, H. J. & White, S. D. M. 1996, MNRAS, 282, 347 
Nagai, D., Kravtsov, A.V., & Kosowsky, A. 2003, ApJ, 587, 524 
Neill, J. D., Hudson, M. J., & Conley, A. 2007, 661, ApJL 123 
Ruhl, Z. ct al., 2004, Proc. SPIE, Vol. 5498, p 11-29 
Sarkar, D., Feldman, H.A. & Watkins, R., 2007, MNRAS, 375, 
69 

Scljak, U., & Zaldarriaga, M., 1996, ApJ, 469, 437 
Stebbins, A., 1997, astro-ph/9705178 

Strauss, M.A., & Willick, J. A., 1995, Phys. Rep. 261, 271 
Sunyaev, R. A. & Zeldovich, I. B., 1980, MNRAS, 190, 413 
Vittorio, N., Juszkiewicz, R., & Davis, M., 1986, Nature, 323, 132 
Wang, L.F. 2007, arXiv:0705.0368 [astro-ph] 
Watkins, R. & Feldman, H.A. 2007, MNRAS, 379, 343-348 
Wright, E., 1998, ApJ, 496, 1 

Yoshida, N., Shcth, R. & Diaferio, A. 2001, MNRAS, 328, 669 
Zcl'dovich, Y.B. & Sunyaev, R., 1969, Ap&SS, 4, 301 
Zhang, P., Stebbins, A., Juszkiewicz, R. & Feldman, H. A., 2004 
arXiv:astro-ph/0410637 



6 REFERENCES 
REFERENCES 

Aghanim, N., Gorski, K.M., & Puget, J. L., 2001, A&A, 374, 1 
Aghanim, N., Hansen, S., & Lagachc, G., 2005, A&A, 439, 901A 
Atrio-Barandcla, F., Kashlinsky, A., & Mucket, J. P., 2004, ApJ, 
601, Llll 

Bardeen, J.M., Bond, J. R., Kaiser, N. & Szalay, A. S., 1986, ApJ, 
304, 15 

Bennett, C.L., et al. 2003, ApJs, 148, 97 
Bhattacharya, S. and Kosowsky, A., 2007, ApJ, 659, L83 
Bonvin, C, Durrcr, R., & Kunz, M. 2006, PRL, 96, 191302 
Courtcau, S., Willick, J., Eds., 2000, Cosmic Flows 1999, ASP 

Conference Series Vol. 201, Sheridan Books, Astron. Soc. Pac, 

San Francisco. 

da Costa, L. N., Nasser, A., Frcudling, W., Giovanclli, R., Haynes, 

M. P., Salzer, J. J., Wegner, G., 1998, MNRAS, 299, 425 
Diaferio, A., et al. 2005, MNRAS, 356, 1477 

Estrada, J., Sefusati, E. & Frieman, J., 2008, arXiv:0801.3485 
[astro-ph] 

Feldman, H., ct al. 2003a, ApJ, 596, L131 

Feldman, H. A., Watkins, R. Melott, A. L. & Chambers, S. W., 

2003b, ApJ, 599, 820 
Feldman, H. A. & Watkins, R., 2008, MNRAS, 387, 825 



© 0000 RAS, MNRAS 000, 000-000 



